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Spatio-spectral properties of twin beams at their transition from low to high intensities are an¬ 
alyzed in parametric and paraxial approximations using the decomposition into paired spatial and 
spectral modes. Intensity auto- and cross-correlation functions are determined and compared in 
the spectral and temporal domains as well as the transverse wave-vector and crystal output planes. 
Whereas the spectral, temporal and transverse wave-vector coherence increases with the increasing 
pump intensity, coherence in the crystal output plane is practically independent on the pump in¬ 
tensity owing to the mode structure in this plane. The corresponding auto- and cross-correlation 
functions approach each other for larger pump intensities. Entanglement dimensionality of a twin 
beam is determined comparing several approaches. 
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I. INTRODUCTION 

The nonlinear process of parametric down-conversion 
[l| is the most frequently used process for the generation 
of light with nonclassical properties. It provides entan¬ 
gled photon pairs in its spontaneous regime 0. Pho¬ 
tons comprising a photon pair can be entangled in differ¬ 
ent degrees of freedom including polarization, frequency, 
emission direction or orbital angular momentum. Entan¬ 
glement in all these degrees of freedom has been found 
useful both for testiM the rules of quantum mechanics 
3 and applications [J. 

On the other hand, parametric down-conversion gener¬ 
ates the so-called twin beams when stimulated emission 
is important. Such twin beams are composed of a signal 
and an idler fields with large intensities that are mutually 
strongly correlated. These correlations occur both in the 
spectra and emission directions as a consequence of the 
properties that give different kinds of entanglement at 
single-photon level. Moreover, the intensity correlations 
are so strong that they violate the standard shot-noise 
limit (sub-shot-noise correlations) @. This nonclassi¬ 
cal property originates in the genuine pairwise character 
of parametric down-conversion at its quantum level. An 
experimental evidence of sub-shot-noise intensity corre¬ 
lations has been given in Such states are useful 

also for quantum imaging. Even imaging based upon 
sub-shot-noise intensity correlations has been recently 
demonstrated 0[ll. 

Spontaneous parametric down-conversion with its pro¬ 
duction of entangled photon pairs has been studied by 
far the most frequently. In theory, the first-order solu¬ 
tion of the appropriate Schrodinger equation provides a 
two-photon amplitude that determines all properties of 
entangled photon pairs [I10. This simple formalism 
is also easily applicable to more complex nonlinear struc¬ 
tures generating photon pairs (waveguides, fibers, lay¬ 
ered structures, Bragg-reflector waveguides, periodically- 
poled structures, for details, see, e.g., [l^). At present, 
photon pairs with more-less arbitrary properties can 


be efficiently generated due to many kinds of available 
sources. 

On the other hand, investigations of twin beams have 
been concentrated to more intense twin beams [T^ be¬ 
cause of the lack of detection techniques capable to de¬ 
tect intensities at the transition from the single-photon 
level to the intense (’classical’) one m. This has re¬ 
quired intense pump lasers that have provided quasi- 
monochromatic beams (usually picosecond pulses). As 
a consequence, the developed theoretical models usually 
invoke the quasi-monochromatic approximation. When 
combined with the pump-field quasi-plane-wave approxi¬ 
mation, a twin beam in the model has been decomposed 
into many more-less independent pairs of signal and idler 
modes localized both in the spectrum and transverse 
wave-vector plane [isl - l^ . The dynamics of individual 
pairs of modes has then been treated by the solution of 
linear Heisenberg equations. Under more general condi¬ 
tions, numerical solution of the Maxwell equations using 
a statistical ensemble of initial conditions has occurred 
useful [^. 


The Schmidt decomposition [23|,l2J| of two-photon am¬ 
plitudes introduced for pure states at single-photon level 
has become popular in the last couple of years due to its 
ability to quantify entanglement in larger Hilbert spaces 
and to reveal the genu ine dual structure of a bipartite en¬ 
tangled state Such modes can even be selected 

from a beam 27H30j or changed on demand. However, 


the revealed paired modes have been found suitable also 
for making the bridge between the perturbation theory 
of weak paired fields and that of the intense beams. The 
corresponding theory has been based upon the solution 
of Heisenberg equations for individual and independent 
paired modes, similarly as the theory developed earlier 
for intense twin beams. Contrary to this theory based on 
localized modes, it uses the Schmidt modes spread over 
the whole spectrum or transverse wave-vector plane [3l[ . 
This makes the theory suitable for describing coherence 
of the twin beams and especially its growth with the in¬ 
creasing pump intensity. It has already been applied to 
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describe spectral properties and amplitude squeezing of 
weaker as well as more intense twin beams |32l-l34| . Also 
ang ular properties of the twin beams have been addressed 
(35l |. Even nonlinear waveguides with back-scattering 
have been investigated using this approach [3^ . 

However, the spectral and spatial properties of a twin 
beam have only been considered separately in this ap¬ 
proach up to now. This represents a serious simplifi¬ 
cation as the spectral and spatial modes are inevitably 
mutually coupled in the nonlinear interaction. The con¬ 
sideration of only the spectral (or spatial) modes allows 
to understand the behavior of twin beams only to certain 
extent as the approach is not able to describe correctly 
a common dynamics of both spectral and spatial degrees 
of freedom. The ’one-dimensional Schmidt-mode models’ 
developed up to now are in fact simple ’fenomenological’ 
models that are conveniently applied for interpreting the 
experimental results obtained under specific conditions 
[strong spatial (spectral) filtering for the spectral (radial 
transverse direction) model]. On the other hand, real 
down-conversion occurs in a nonlinear crystal with all 
possible spatio-spectral modes participating in the inter¬ 
action. As bulk nonlinear crystals are commonly used 
for the generation of intense twin beams, the number 
of participating modes is large. Moreover, the number 
of modes giving an important (intensity) contribution to 
the generated twin beam is also large. That is why, a 
general spatio-spectral model of twin-beam generation is 
necessary. 

In this contribution, we develop the Schmidt-mode 
approach for describing such a general spatio-spectral 
twin beam. Using the paired spatio-spectral modes, we 
address coherence of the twin beams determining their 
auto- and cross-correlation functions in the spectral and 
temporal domains, transverse wave-vector plane, and 
crystal output plane. Properties of the twin beams man¬ 
ifested under different experimental conditions are com¬ 
pared. Special attention is paid to the dependence of 
coherence on the pump intensity. Coherence properties 
of the twin beams together with their photon-number 
statistics are used to compare several suitable quantifiers 
of their dimensionality [371138|| . This provides a complete 
picture of an intense spatio-spectral twin beam. 

The developed model can easily be generalized to more 
complex nonlinear structures with nonlinearity homoge¬ 
neous in the transverse plane including poled nonlinear 
crystals [H, and nonlinear layered structures [l^ . It 
can even be applied for the description of intense twin 
beams originating in the nonlinear process of four-wave 
mixing, both in nonlinear crystals with nonlinearity 
and atomic ensembles. The four-wave mixing in atomic 
ensembles [4l|, though being partly noisy [d^, is very 
attractive due to the high effective nonlinear coup ling 
constants. Noiseless spatially-resolved amplification 
as well as entangled images [Q in this scheme using 
vapors have been experimentally demonstrated. 

Recently, the first experimental investigations of ultra¬ 
intense twin beams have been reported both for multi¬ 


mode [13, Id^ - ldSj j as well as single-mode (bright squeezed- 
vacuum states) fields [d^. The effects of pump-field de¬ 
pletion have been observed. Also back-flow of energy 
from the twin beam into the pump field may occur. These 
processes affect spectra as well as transverse profiles of 
the twin beams [d5l - ld7j | . Here, we restrict our attention to 
the case of un-depleted pump fields (parametric approxi¬ 
mation). However and importantly, the generalization of 
the theory to account for pump-field depletion and back- 
flow of energy is possible due to the fact that the model 
incorporates all spatio-spectral degrees of freedom. This 
generalization will be reported elsewhere. 

The paper is organized as follows. A spatio-spectral 
model of parametric down-conversion based upon the 
Schmidt paired modes is developed in Sec. H. Quanti¬ 
ties characterizing spectral and temporal properties of 
twin beams are defined in Sec. HI. Several suitable quan¬ 
tifiers of dimensionality of a twin beam are introduced 
in Sec. IV. Spectral and temporal properties of twin 
beams as functions of the pump intensity are discussed in 
Sec. V. Properties of twin beams in the transverse wave- 
vector plane and the crystal output plane are analyzed 
in Sec. VI. Conclusions are drawn in Sec. VH. In Ap¬ 
pendix A, twin beams are described in their transverse 
wave-vector plane (far field) and the crystal output plane 
(near field). 


II. THEORY OF A SPATIO-SPECTRAL TWIN 
BEAM 


Optical parametric down-conversion and its evolution 
along a nonlinear medium characterized by tensor d of 
second-order nonlinear coefficients is described by the 
momentum operator Gjnt written in the interaction rep¬ 
resentation as follows (Mini: 


G'int(z) = 4eo / dxdy 


dt 


[d : 4+) (r, t)4(-) (r, t)(r, t) + H.c.] ; (1) 


r = {x,y,z). In Eq. ([T]), describes the positive- 

frequency part of a classical pump electric-field amplitude 
and Es~^ stands for the negative-frequency part 

of a signal- [idler-] field operator amplitude. Symbol : is 
shorthand for tensor shortening with respect to its three 
indices, cq is permittivity of vacuum and H.c. replaces 
the Hermitian conjugated term. 

The electric-field amplitudes [ui''’^ (r, t) = 

of the interacting fields can be decomposed 
into harmonic plane waves with wave vectors and fre¬ 
quencies UJa'- 


£;(+)(r, t) = 3 [ dka El+\ka) exp(ikar - iujJ), 


a = p, s, i. 


( 2 ) 
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We assume that the interacting fields can be described 
in paraxial approximation which provides the relation 
k = ikx,ky,k^) = {kx,ky,k - [kl + k‘^]/2k), k = 

kx + ky + kl, valid for fields propagating close to the 
z direction. 

We consider a strong pump field with the Gaussian 
spectrum and Gaussian transverse profile. It is described 
in paraxial approximation as follows: 


H is the reduced Planck constant. We note that Ua = 
cka/uJa gives the index of refraction of field a. The cre¬ 
ation and annihilation operators fulfil the usual boson 
commutation relations appropriate for the quantization 
of photon flux [H, , 

[da(k^,Wa),a|;,(k';|i,w;,)] = (5aa'<5(k^ - k';[;) 

X (5(a;a - Wa'), (9) 


4+)(r,t) = 


1 








X e'xp{ikp^xx) e'xp(ikp,yy) exp(zfcpz) 
X exp —i 


k"! „ k'i ,, \ 

— z ) exp(-ia;pt); (3) 




kp = (kp^xjkp^y) and fcp = fcp(a;p). Temporal spectrum 
£'||(a;p) of the Gaussian pump pulse is expressed as: 






(4) 


where S means the Dirac 6 function and Saa' stands for 
the Kronecker symbol. 

Substituting expressions ([3]) and (O into Eq. ([T]) for 
momentum operator Gint we arrive at the formula 

Gint{z) = f dk^ f dk^ rdw, HdoJi 

J J JO Jo 

poo 

dujp 6{ujp -uJs- uJi)El{u}p)^^TL{k^,kl~) 

X exp (i[fcp(a;s -|- Wi) - ^^(ws) - fci(a;i)]z) 

X a^(k^,a;s,z)a|'(k/-,a;i,z) -hH.c., (10) 


According to Eq. (|T|), the pump pulse has amplitude ^p, 
duration Tp and carrying frequency Wp. Provided that 
the pulsed pump field source has power Pp and repetition 
rate /, amplitude ^p is given by the relation 


^P = 



(5) 


where c is the speed of light in vacuum. The pump field 
radially symmetric in its transverse plane is characterized 
by the Gaussian spatial spectrum 


£;p^(k^) = ^exp 


Wlik^x + k^y) 


( 6 ) 


where Wp gives the beam radius. 

The signal and idler electric-field operator amplitudes 
Es ^ and E^ ^ are decomposed similarly as the pump 
field in paraxial approximation: 


[ dki r dWaEi-\ki,CUa) 
j do 

X exp(-ffca,£ca:) exTp{-ika,yy) exp(-ifcaz) 

/ ^2 _|_ j^2 \ 

X exp exp{iu}at), a = s,i. (7) 


The spectral operator amplitudes Ea ^(k^,Wa) can be 
expressed using creation operators d^(k^,a;a) that add 
a photon into the mode with transverse wave vector k^ 
and frequency uJa- 


Ei ^(k^ ,a;„) = -G 


2eoc^ka 


ai(kf,a;a); (8) 


where deS is an effective nonlinear coefficient. Eunction 
Tl describes correlations between the signal and idler 
fields in the transverse wave-vector plane: 


7’L(k 


S , k; ) = / dkp 5(kp - k, - k; )£;p (kp ) 


X — y dz exp I —i 


.J-|2 


2kn 


,J-|2 


2ks 


Ikf 


±12 


2h 


;(ii) 


|ko P = ^a.xd-k'^ y. We note that an average value of the 
phase mismatch determined along the crystal of length 
L occurs in formula (HU. 

These correlations are conveniently expressed using 
dual ortho-normal transverse modes of the signal and 
idler fields. These modes are revealed by the Schmidt de¬ 
composition of the normalized function T^, Tl = t-^T^ 
and = J dk^ f c?k[*-|TL(k^, k[*-)p. As we are in¬ 
terested in the radially symmetric geometry, the use 
of radial variables k^ and (pa is convenient [k^ = 
(/c^ cos((pa), fca sin((/ 2 Q))]. We note that the considered 
radial symmetry is broken for narrow pump beams owing 
to the crystal anisotropy . In the radially symmet¬ 

ric geometry, the normalized function T£ can be rewrit¬ 
ten into the form: 

^ OO 

EE(k^,k^) = ^ TL,m(kE,kJ)explim(<ps-¥’i)]- 

m— — oo 

( 12 ) 

Eunctions TL^m introduced in Eq. m and defined as 

p2‘K 

TL,m{kE,kJ) = / d(v5s - <Pi)TE{kE,kE) 

Jo 

X exp [-im{(ps - v^i)] (13) 
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can be decomposed as follows: 

- “ 

Y' k^k^TL,m{ki, k^) = \^iUs^rnl{ki)Ui.,nl{k^). 

1=0 

(14) 

Eigenfunctions Us,mi and Ui^mi form the ortho-normal 
dual bases and denote the corresponding eigenval¬ 
ues. Substituting Eq. ([T4)l into Eq. (lED we reveal the 
Schmidt decomposition of the normalized function T£: 


OO OO 


\Jk^k^T£{ki-,ip^, k^, ifi) = Y 


m——oo l—O 
u-L 


^ ti,ml (^i 5 ■ 


(15) 


The transverse mode functions ts,mZ and occurring 
in Eq. dTK]) take the form: 


is,ml {kg , IPs') — 


ii,ml {k\ , p^i) — 


Us,mi(k^)exp(imips) 
Ui^mi(k[‘-) exp(-im(pi) 


(16) 


Introduction of field operators aa,mi(oJa, z) related to 
transverse mode functions ta,mi, 

poo n 27 T 

O‘a,ml(0JaT z) = I dk^^ / dipata^mlO^a t ^a) 

Jo Jo 

X daik^,(Pa,UJa,z), O = S, i, (17) 

allows to rewrite the interaction momentum operator 
Gint in Eq. m as follows: 


_ 2Sdefft'‘‘ .^_L f°° . 

Gint{z) — - 3 / . ^ml I 

c2 „ , Jo 


duJi 


Vhh 


X eI|(ws -I- Wi) exp (i[fcp(a;s -I- Wi) - fcs(ws) - ki{uj\)\z) 

X z)dYi{u:p z) + H.c. (18) 

If the nonlinear interaction is weak, we can obtain a 
perturbation solution of the corresponding Schrddinger 
equation and express the output state |'0)out in the form: 


dz Gint('^)|V^)im 


(19) 


where |'!/')in is the input signal and idler state. Substi¬ 
tution of Eq. (IT^ into Eq. (IT^ and consideration of the 
input vacuum state |vac) result in the formula 


m,l 


■ml 


dujs / duJiFL{uJs,uJi) 


X (‘*'1:0)1 vac), 


( 20 ) 


where 


X / dz exp {i[kp{uis + uji) — ks{u!s) — ki{uji)]z). (21) 
Jo 


We note that the vacuum state |vac) is omitted in the 
expression for the output state |t/’)out in Eq. (EOl). 

Using eigenfunctions fs,q and fi^q and eigenvalues a| 
of the Schmidt decomposition of the normalized function 
Pl [Pl = = f dujsf duJi\FLiu}s,0Ji)\'^], we can 

rewrite Eq. m as follows: 


FL(a;s,Wi) = /II^A|/s,,(ws)/i,q(a;i). (22) 

9=0 

New field operators da,miq defined as 

poo 

0^a,mlq — I duJa f 0); Q- = S, i, 

(23) 


provide a simple formula for the output state |^)out: 

|V')out = Y ^^l^l^lmlq^lmlq\^^(^)- (24) 

m,l,q 

According to Eq. (IMl) the output state |V’)out is com¬ 
posed of photon pairs in independent paired modes num¬ 
bered by indices {m,l,q) with probability amplitudes 

Using the paired modes revealed both in the transverse 
wave-vector plane and spectrum we rewrite the ’aver¬ 
aged’ momentum operator dzGintiz)/L from Eq. (fTSl) 
into the form: 


Gt:dz)=-'-^ Y fl >^ml>^l<mlqi^)^lmlqi^) 

m— — oo l,q—0 

-l-H.c. (25) 

using the operators da,miq defined in Eq. (j23p . The cru¬ 
cial advantage of the ’averaged’ momentum operator Gfj((. 
is that it ’diagonalizes’ the nonlinear interaction leav¬ 
ing the separated Heisenberg equations for each pair of 
modes. We note that some of the paired modes are degen¬ 
erate in certain symmetric configurations (e.g. collinear 
spectrally-degenerate emission) in the sense that both 
the signal and idler photons are emitted into the same 
spatio-spectral mode [H^. However, we do not consider 
explicitly such modes here. Considering an {m,l,q)-th 
mode, the Heisenberg equations are written as follows: 

(26) 

using effective nonlinear coupling constants Kmiq, 


fC _ y-L II 

i^mlq — j -^ml-^q- 


L 


( 27 ) 
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The solution of linear equations (l26)) for an (m, I, (i')-th 
mode and the crystal of length L takes a simple form: 

ds,m/g(T) = COS\l(^KjnlqL^^s,mlq(S^) 

( 0 ) 

+ smh{KmigL)al^igiO). (28) 

The solution for a given transverse mode (m, 1) can be 
conveniently expressed in the matrix form: 


where Ska gives the declination of the transverse wave 
vector of field a. The derived formulas valid for the close- 
to-collinear geometry can be applied in general also in 
the non-collinear case provided that the following formal 
substitution is used: 

ka <— ka cos(d°), Ska <— Ska cos(d°)^. (34) 

III. SPECTRAL AND TEMPORAL 
PROPERTIES OF TWIN BEAMS 


^s,mz(T) — (0) “t” ^ 

a-i,mz(T) = Ui^mzai_mi(0) + ,^j(0). (29) 

The matrices Us,mi, Ui,mi and Vmz introduced in 
Eq. p9l) are written in their singular-valued decompo¬ 
sitions as follows: 


Ua,mi — ^a,ml-^ml^a,mh ^ h 

\r _ IP \V ipt 

^ ml — ^ s,ml-^^rnl^ i^ml' 


(30) 


Columns of the matrices Fs,mi (Fi,mi) in Eq. (1501) are 
given by eigenmodes fs,q ifi,q) of the Schmidt decompo¬ 
sition written in Eq. Elements of the diagonal ma¬ 

trices A^; and A))); are derived from the solution given 
in Eq. (1^ , 

^ml,qq — ^mlq — COSh(AmigT), 

^ml^qq — ^mlq — sinh(AmigT). (i^^) 

The solution (l28l) allows to derive the mean values of 
experimental physical quantities. Spectral and tempo¬ 
ral quantities are determined in Sec. 3 below. Spatial 
quantities are defined in Appendix A. Numbers of modes 
constituting the twin beam are described in Sec. 4. 

We note that the numerical results obtained in [s^ 
show that mild broadening of the modes determined from 
the perturbation solution of the Schrodinger equation oc¬ 
curs for strong pumping of the nonlinear process. 

We also note that, in the considered radially symmet¬ 
ric non-collinear geometry with the pump field at normal 
incidence, the signal and idler fields propagate along the 
radial emission angles ds and di, respectively. The cen¬ 
tral radial emission angles d® and d? corresponding to 
the central frequencies w® and urf are given by the con¬ 
servation of energy and transverse wave vectors: 

= a;° - sin(d°) = k° sin«), (32) 

= fca(a;°). The central transverse wave vectors 
are then given as cos(i5°), a = s, i. Paraxial ap¬ 

proximation along the radial emission angle provides 
the following formula for wave vector (a = s, i): 


k„ = 


([■ 


k^° + Ska] cos{(pa), [fca ° + Ska] sin(i^o). 
Ski cos(ds)^ 


ka - 


2ka 


cos(d°)), 


(33) 


We assume that the transverse profiles of twin beams 
are not experimentally resolved and so the experimental 
mean values are obtained by averaging over the trans¬ 
verse modes. Then the signal-field intensity spectrum 
ris^cj is expressed as follows: 

ns,ui(^s) = (a|[(a;s,L)as(a;s,L))± 

= (35) 

ml q 

Symbol {)± denotes quantum mechanical averaging com¬ 
bined with averaging in the transverse plane. The num¬ 
ber Ns of generated signal photons is determined along 
the formula 

pCO _ 

Ns = duJs ns,u,{uJs) = E E (36) 

ml q 

Averaged signal-field intensity spectral correlations are 
characterized by the fourth-order correlation function 
given as: 

As,a; (^s, a^s) — T)] 

X A[a;[(w(,L)as(a;(,L)] :)_l 

= Y.\AlmiA^s.uj'^)\"- (37) 

ml 

The signal-field amplitude correlation function ^ 

belonging to mode (m, 1) is written in the form: 

As,ml,uji^ST^s) ~ (a|(Ws, T)as(Ws, T))_L,mi 

= Y.flq{^.)U,qU)Vlig. (38) 

<? 

Intensity spectral cross-correlations between the signal 
and idler fields are quantified by the following fourth- 
order correlation function: 


Cu (a^s, ^i) 


(A/": A[ai(a;s,T)as(ws,T)] 

X A[a((wi,L)ai(wi,L)] :)_l 

2 


E 


^ ( fs,q{^s] fi,qi^i)UijilqNijilq 
Q 


(39) 


Temporal electric-field amplitude and intensity corre¬ 
lations in the twin beams outside the nonlinear crystal 
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can be expressed, similarly as their spectral correlations, 
in terms of temporal eigenfunctions fa,q{ta) determined 
by the Fourier transform: 



The averaged signal-field photon flux is then derived 
in terms of functions fs,q, 

Is,tits) = 2 eoc{E^~\r^,L,ts)Ei'^\r^,L,ts))± 

= (41) 

ml q 

The averaged signal-field intensity temporal correlation 
function is expressed similarly as the spectral corre¬ 
lation function As,uj given in Eq. (EH), 

As,tits,t's) = (2eoc)^(A/': (r^, L, 4) 

X Ei+\r^,LAs)]A[&-\r^,LX)^+\r^,LX)] :)± 

^'EKrni,t{tsXs)\^ ■ (42) 

ml 

The signal-field amplitude temporal correlation function 
A^mi t mode (m, 1) is determined along the formula 

Al^iXt^,t's) = 2eoc{^-HrXL,ts)Ei-^HrXL,ts))±,mi 

= E/s%(4)/.,,(i()14^i,- (43) 


Also the averaged intensity temporal cross-correlations 
between the signal and idler fields can be quantified in 
the same vein as in Eq. (15^ : 


CtitsAi) = ( 2 eoc)"(AA : A[Et\rXL,ts) 

X EXHrXL,ts)]A[Mx\rXL,ti)M^\rXL,ti)] :)± 


= E 


ml 


^ ^ fs,q(j^s')fijqiiijUmlq^mlq 


(44) 


As the number of transverse modes is usually large, 
their eigenvalues XX form quasi-continuum. In this case, 
we may introduce a suitable probability function g\ and 
make the following replacement in the above formulas: 


E 



dX-^gx{X-^). 


(45) 


This makes the numerical computations considerably 
faster. 

The mode {m,l,q) = (0,0,0) having the largest value 
of the product of the Schmidt eigenvalues becomes 
dominant in the limit of large pump power (Pp —> oo) in 
the used un-depleted pump approximation. The spectral 
characteristics ns,u;, As,uj and Cs,u) then attain the simple 
form: 

^s,a;(^s) — A^s I/*s,0 (^s) I ? 

As,Xu^s.X) = N^Afs,,Xs)?\fs,,{X)\\ 

(^(.^(Ws, Wi) = Aqq UgQ QlQg_Q|/s,o(Ws)| |/i,o(Wi)| , 

( 46 ) 


where Ns = g gives the number of emitted signal pho¬ 
tons. According to Eqs. (HSl) . the twin beam is spectrally 
composed of independent single-mode signal and idler 
fields in this high-intensity (classical) limit. We note that 
one dominant paired mode constitutes the twin beam also 
in the transverse wave-vector plane and the crystal out¬ 
put plane. So the signal and idler fields are spatially and 
spectrally independent but internally fully spatially and 
spectrally coherent. 

Similar quantities as defined above in the spectral and 
temporal domains are used for describing the twin beams 
in their transverse wave-vector plane and crystal out¬ 
put plane. Modes in the crystal output plane are de¬ 
termined from those of the transverse wave-vector plane 
using the Fourier transform, similarly as the temporal 
modes have been derived from the spectral modes. The 
radial symmetry of twin beams results in harmonic az¬ 
imuthal modes in the crystal output plane. It also pro¬ 
vides the radial modes in the crystal output plane de¬ 
termined from those of the transverse wave-vector plane 
using the transformation based on the Bessel functions 
(for details, see Appendix A). 


IV. DIMENSIONALITY OF THE TWIN BEAM 


Dimensionality of a twin beam can be determined ei¬ 
ther using its paired properties or properties of the indi¬ 
vidual signal and idler fields. In the first case, dimension¬ 
ality of entanglement is obtained. In the second case, the 
numbers of independent modes constituting the signal (or 
idler) field and defined in statistical optics are reached. 
Entanglement dimensionality for a general noisy twin 
beam is quantified via negativity [s^. Considering pure 
states of the noiseless twin beams, the Schmidt number 
can be applied f or q uantifying entanglement dimension¬ 
ality as well [H, l24j . This number can even be reached 
without making the Schmidt decomposition The 

general formulas can be recast into a simple form for 
quasi-monochromatic or quasi-homogeneous fields M- 
Compared to weak fields, the analysis of intense twin 
beams is more difficult, as the decompositions not only 
in the spatial and spectral domains but also in the 
Hilbert spaces of individual paired spatio-spectral modes 
spanned by the Fock-number states would be needed. 
That is why, we apply here a simpler approach for deter¬ 
mining entanglement dimensionality based upon defining 
creation operators for photon pairs (for details, see [3^1. 
Entanglement dimensionality K of the twin beam is ob¬ 
tained in this approach as follows: 


K = 


L L^mlq ^ mlq ^ mlq 
Z^mlq mlq mlq 


(47) 


We note that formula (|T7)) reduces to the usually used 
Schmidt number of spatio-spectral modes for weak twin 
beams. 
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Formula (gzl) can also be applied to provide aver¬ 
age number of effectively populated paired spectral 
modes: 

= ^ m VA ■ (48) 

ml Z^q'^mlq'^mlq 


In Eq. gzl), Pmi gives the probability of having a photon 
pair in mode (m, 1): 


y v^, 

_L mlq 

Pml — ■ 

L^mlq ^ mlq 


(49) 


Similarly as in the spectrum, average number Kkip of 
effectively populated modes in the transverse wave-vector 
plane is obtained along the formula: 


(s 


r/^ 1/2 

E ll \^-^ml mlqmlq ^ 

Pq fjZ 773 

g mlq ^mlq 


(50) 


ns,mlq — ^\,mlq^^)^s,mlq{L) ■ 

Also, the formula for averaged number AT” of spectral 
modes can be written, in analogy with the derivation of 
Eq. (H51) from Eq. (ITTl) : 




(53) 


_ 


(S 


q y^s,mlq 


( (A^ : ^^s,mlq '■) {ns,mlqY^ 


yqVLq 


The averaged number of modes in the transverse 
wave-vector plane can be determined by a formula anal¬ 
ogous to that written in Eq. (l5^ [compare Eqs. (l48l) and 

gni- 

The ratio of the width Arig of a signal-field inten¬ 
sity profile and the width AA^ of the appropriate signal- 
field amplitude auto-correlation function. 


using probability of having a photon pair in mode q\ 

y 

pI = . (51) 


'T ■ 

y^mlq mlq 


The numbers K^j and Kk^p of paired modes in the spec¬ 
trum and transverse wave-vector plane, respectively, can 
alternatively be determined as the ratio of the width An^ 
of, say, the signal-field intensity profile and the width 
of intensity cross-correlation function AC in the appro¬ 
priate variable. This ratio known as the Fedorov ratio 
|62l | coincides with the number AT^ of paired modes given 
in Eq. (1481) for weak twin beams with a Gaussian two- 
photon amplitude [^. It has been shown in that 
both numbers are close to each other for general weak 
twin beams. 

Modes in the signal and idler fields are ideally paired as 
well as the signal and idler photons in individual spatio- 
spectral modes for the considered noiseless twin beams. 
That is why, dimensionality of the twin beam can also be 
determined from the number of modes and their popula¬ 
tions counted either in the signal or idler field. Applying 
the coherence theory [b^ , an effective number of indepen¬ 
dent modes (degrees of freedom) in the signal (or idler) 
field can be obtained from its photon-number statistics. 
The resulting number AT" of modes constituting, e.g., the 
signal field is given by the formula valid for a multimode 
thermal field |66| : 


AT” = 


(ymlq('^s,mlq}^ 


yiimlq ((A^ ■ ^1,mlq ') {ns,mlqY^ 

E. 


(y E^ V 

\ A^mlq ^ mlq J 
y^mlq * mlq 


(52) 




AA^ 


(54) 


defined in any variable represents also a good quanti¬ 
fier of the number of independent modes of a twin beam 
in this variable. We compare different quantifiers of di¬ 
mensionality of the twin beam under real experimental 
conditions below. 


V. SPECTRAL AND TEMPORAL PROPERTIES 
OF INTENSE TWIN BEAMS 

In the numerical analysis, we consider a BBO crys¬ 
tal 8-mm long cut for non-collinear type I process (eoo) 
for the spectrally-degenerate interaction pumped by the 
pulse at wavelength Ap = 349 nm with spectral width 
AAp = 0.1 nm, transverse profile with radius iCp = 1 mm 
and repetition rate / = 400 s“^. This pulse is provided 
by the third harmonics of the Nd:YLF laser at wave¬ 
length 1.047 /rm. Assuming the pump field at normal 
incidence, the signal and idler fields at the central wave¬ 
lengths Ag = A? = 698 nm (i9bbo = 36.3 deg) propa¬ 
gate outside the crystal under the radial emission angles 
■dg = i?? = 8.45 deg. As this configuration is symmetric 
for the signal and idler fields, we further discuss only the 
properties of signal field. We assume that the conditions 
are such that the spectral and spatial properties of the 
twin beam factorize. 

The generated twin beam is composed of roughly 80 
thousand transverse modes at low intensity. It contains 
34 modes in radial direction and 2350 modes in azimuthal 
direction (for more details, see 0). As the number 
of transverse modes is large, we can introduce quasi¬ 
continuum of the Schmidt eigenvalues with its proba¬ 
bility function g\ defined in Eq. (H51) . The probability 













FIG. 1. (a) Probability function g\ of eigenvalues in the 
transverse wave-vector plane and (b) spectral eigenvalues A,; 
trip = 1 X 10“® m, AAp = 1 X m. 



Pp (10'® W) 


FIG. 2. Number Ns of emitted signal photons (solid curve 
with A), number Ns,k'^ of emitted signal photons per mode 
defined by photon-number statistics (solid curve with *), 
number Ns k of emitted signal photons per mode given by 
Eq. (1541) (plain curve), and gain g (dashed curve) as func¬ 
tions of pump power Pp; log denotes the decimal logarithm; 
Wp = 1 X 10“® m, AAp = 1 X 10“^° m. 


function q\ is plotted in Fig. [IJa). It reflects the fact 
that the smaller the eigenvalue the larger the number of 
such eigenvalues. There occur around 80 independent 
spectral modes in the low-intensity regime, as shown in 
Fig- mb). 

The number Ns of emitted signal photons increases 
roughly exponentially with the increasing pump power 
Pp for more intense fields [l^, as shown in Fig. [21 The 
curves in Fig. |2] giving the number of emitted signal pho¬ 
tons per one mode show that the probabilities of spon¬ 
taneous and stimulated emissions of a signal photon (to¬ 
gether with its idler twin) become comparable for pump 
powers Pp around 0.5 mW. Exponential increase of the 
number Ns of emitted signal photons occurs already for 
pump powers Pp one order of magnitude lower. It is use¬ 
ful to define gain g of the nonlinear interaction and to 
use it instead of the pump power Pp when comparing 
the theoretical quantities with their experimental coun¬ 
terparts. The gain g arises from a simplified model of 
the nonlinear interaction that assumes only one effective 
mode and the initial vacuum state of the signal and idler 
fields. Formulas (l?7)) . (05)) and (1551) of Sec. II provide in 



FIG. 3. Spectral entanglement dimensionality (plain solid 
curve), number of modes determined from photon-number 
statistics (solid curve with >i=) and number of 

spectral [temporal] modes given by Eq. (1541) (dashed curve 
[dashed curve with A]) as functions of pump power Pp; 
Wp = 1 X 10“® m, AAp = 1 X 10“^° m. 

this case the expression {Xqq = Aq = 1) 

fV, = sinh(t^/ll/fF)2_ (55) 

It suggests the following approximative formula for fitting 
the experimental dependence of the number Ns of signal 
photons: 

Ns = NsfiSinh{g)^, (56) 

where g = go\/P^ and Pg.o and go are suitable con¬ 

stants. The values of gain g assigned to pump powers Pp 
are plotted in Fig. [21 They show the advantage of this 
parametrization: Stimulated emission of photon pairs be¬ 
gins to dominate over spontaneous emission for the values 
of gain g around one and the transition from quantum 
to classical regimes (mesoscopic regime) occurs for the 
values of g around 10. 

Spectral entanglement dimensionality Ki^ determined 
by formula (i48l) decreases with the increasing values of 
pump power Pp (see Fig. Ej). The number of 
spectral signal-field modes as well as the number of 
temporal signal-field modes given in Eq. ([SI by the ra¬ 
tios of appropriate widths and plotted in Fig. [3| are lower 
than the spectral entanglement dimensionality Ki^. The 
comparison of curves in Fig.Ejshows that the experimen¬ 
tally available values of and can successfully be 
used for quantifying dimensionality of the twin beam, to¬ 
gether with the theoretical entanglement dimensionality 
Ki^. The number of modes constituting the twin beam 
can also be derived from the photon-number statistics in 
the signal (or idler) field [lo, 113, IHl ■ In this case, the 
number of modes is given by formula (1531) . It provides 
systematically greater numbers of modes, as the curves 
in Fig. E] show. The values of entanglement dimensional¬ 
ity and number PT" of modes nearly coincide in the 
low-intensity regime. This immediately follows from the 
comparison of Eqs. (HSI) and (1551) in the limit Umiq ~ 1. 
Also, the nnmbers and Ff" of modes equal to 
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FIG. 4. (a) Width Aus^^j of signal-field intensity spectrum 
as a function of pump power Pp and (b) spectrum for 
Pp = 1 X 10“^ W (plain curve) and Pp = 2 x 10“^ W (solid 
curve with *); Wp = lxl0~® m, AAp = 1x10”^'^ m. Spectrum 
ns,u; is normalized according to f diUsns,aj((j-’a)/(j-’s = 1- 


FIG. 6. (a) Width A/s.t (FWHM) of signal-field photon flux 
as a function of pump power Pp and (b) photon flux 7s,t for 
Pp = 1 X 10“^ W (plain curve) and Pp = 2 x 10“^ W (solid 
curve with *); nip = 1 x 10~® m, AAp = 1 x 10”^'^ m. The 
curves in (b) are normalized such that f dts7s,t(is) = 1. 



(a) (b) 

FIG. 5. (a) Widths A^s,i.j of signal-field intensity auto¬ 

correlation function (solid curve) and ACs,!^ of intensity 
cross-correlation function (dashed curve) as functions of 
pump power Pp. In (b), intensity auto-correlation function 
4ls,tj(‘^s) = Als,i.u(iUs,oJs)Ma,u(n)s,oJs) and cross-correlation 
function Cs^i^(uJs) = Cs,uj(oJs,‘^°)/Cs,uj(oJb,‘^?) are plotted for 
Pp = 1 X 10“^ W (plain curves) and Pp = 2 x 10“^ W 
(nearly coinciding curves with >i=); Wp = 1 x 10“® m, AAp = 
1 X 10”“ m. 


the entanglement dimensionality 77^ in the high-intensity 
limit (Pp — >■ cxd). This occurs because the strongest mode 
completely dominates over the other modes in this limit. 

Decrease of the number of signal-field modes with 
the increasing pump power Pp originates in the behav¬ 
ior of spectral widths Aus^uj and Whereas the 

spectral width Arig,^ of the signal-field intensity pro¬ 
file decreases with the increasing pump power Pp [see 
Fig.SKa)], the width AAg of signal-field amplitude auto¬ 
correlation function increases [for the width AAs,tj of in¬ 
tensity auto-correlation function, see Fig. [SKa)]. This 
occurs because the spectral modes with greater eigen¬ 
values Xq become more and more important with the 
increasing pump power Pp. Hand in hand, the role of 
modes with small eigenvalues Xq is suppressed. As the 
modes with large eigenvalues Xq are localized more in the 
middle of the spectrum (for more details, see, e.g., [^). 
narrowing of the signal-field intensity spectrum is natu¬ 
rally observed. This is accompanied by reshaping of the 
spectrum that looses small oscillating tails present 


in the low-intensity regime [see Fig.|4l[b)]. As the genera¬ 
tion of photons by stimulated emission increases with the 
increasing pump power Pp, coherence in the twin beam 
increases. This leads to broadening of the widths AAs,cj 
and AC's,,, of intensity auto- and cross-correlation func¬ 
tions [ 43 . As the curves in Fig.[Sl[b) drawn for two differ¬ 
ent pump powers Pp show, the intensity auto-correlation 
function Ag^ui is wider than its cross-correlation coun¬ 
terpart Cg^uj in the low-intensity regime (for the expla¬ 
nation, see [i^). When stimulated emission begins to 
dominate over spontaneous emission [see Fig. |5l[a)], the 
auto-correlation and cross-correlation functions approach 
each other. In the high-intensity limit Pp —>■ 00 , the sig¬ 
nal and idler fields are single-mode and so they are spec¬ 
trally coherent. 

The signal field is emitted in the form of a short 
pulse. It is composed of the temporal modes /s,g given 
in Eq. (gni)- These modes behave similarly as their spec¬ 
tral counterparts [s^. Thus, a g-th mode has q max¬ 
ima and q — 1 zeroes in its intensity temporal profile. 
Also, the greater the number q, the wider the mode. 
The field transition to the high-intensity regime looks 
as follows. The signal pulse is in general l ong er than the 
pump pulse in the low-intensity regime 1^. However, 
as shown in Fig. EKa) the signal pulse shortens with the 
increasing pump power Pp due to the nonlinear inter¬ 
action described in momentum operator Gjnt written in 
Eq. II]). The signal pulse is also delayed with respect 
to the pump pulse [see Fig. IHb)] as a consequence of 
different group velocities of two pulses inside the crystal 
[T^ . Coherence in the signal field as well as coherence be¬ 
tween the signal and idler helds increase with the increas¬ 
ing pump power Pp due to stimulated emission, as docu¬ 
mented in Fig.|7][a). The intensity auto-correlation func¬ 
tion Ag^t is narrower than the intensity cross-correlation 
function Cg^t in the time domain and low-intensity regime 
[see Fig. l2l[b)]. This is opposed to the behavior of spec¬ 
tral correlation functions. It originates in properties of 
the Fourier transform. The cross-correlation and auto¬ 
correlation functions are close to each other for greater 
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FIG. 7. (a) Widths AAs,t (FWHM) of signal-field intensity 
auto-correlation function (solid curve) and ACa,t of inten¬ 
sity cross-correlation function (dashed curve) depending on 
pump power Pp. In (b), intensity auto-correlation function 
= As.t(ts,ir“) and cross-correlation 
function C's_t(ts) = C's,t(ts, are plotted 

for Pp = 1 X 10“’’ W (plain curves) and Pp = 2 x 10“^ W 
(nearly coinciding curves with *); and 1™“ give the times 
with maximal photon fluxes in the signal and idler fields, re¬ 
spectively; Wp = 1 X 10“® m, AAp = 1 X 10“’° m. 

values of pump power Pp, as shown in Fig. EKb). In the 
high-intensity limit Pp —>■ oo, the twin beam is found in 
a separable temporally coherent state composed of the 
signal- and idler-field temporal modes /s,o and fi^ writ¬ 
ten in Eq. (PO]). 


VI. PROPERTIES OF INTENSE TWIN BEAMS 

IN THE TRANSVERSE WAVE-VECTOR AND 
CRYSTAL OUTPUT PLANES 

We analyze properties of the twin beam in the wave- 
vector transverse plane (far field) and the crystal output 
plane (near field) assuming spectral (or temporal) aver¬ 
aging. 

Entanglement dimensionality in the transverse 
plane gives around 80 thousand modes in the low- 
intensity regime. It decreases with the increasing pump 
power Pp (see Fig. |S]) [13, This behavior is similar 
to that found in the spectral and temporal domains. It 
can be explained in the same way. Entanglement dimen¬ 
sionality Kkip in the transverse plane and number 
of signal-field transverse modes provided by the photon- 
number statistics are close to each other, as shown in 
Fig. [8l These numbers can alternatively be experimen¬ 
tally estimated using the product of ratios of in¬ 

tensity widths and widths of amplitude auto-correlation 
functions both in radial and azimuthal transverse wave- 
vector directions applying formula (j54[) . Factorization of 
the number of modes into its radial and azimuthal con¬ 
tributions is valid in our geometry in which the photons 
are emitted into a narrow ring in the wave-vector trans¬ 
verse plane. In the low-intensity limit, there is around 34 
[2350] modes in radial [azimuthal] wave-vector direction. 
Around 10 [1000] modes are found in radial [azimuthal] 
wave-vector direction for the pump power Pp = 50 mW. 



FIG. 8. Entanglement dimensionality Kkp (plain solid 
curve), number of signal-field modes determined from 
photon-nnmber statistics (solid cnrve with *) and num¬ 
ber K^).^ of signal-field modes in the transverse 

wave-vector [crystal ontpnt] plane given by Eq. (|54l) (plain 
dashed curve [dashed curve with A]) as they depend on 
pump power Pp. Number equals K^r given in 

Eq. (1541) in which the width AA® j. = 2 dva (tb — 
’'0^s,r-(’'s,r()//)rdraAJ^^(rs,r(); Wp = 1 x 10“° m, AAp = 
1 X 10“’° m. 


On the other hand, the signal and idler photons form a 
disc centered around x = y = 0 m in the crystal out¬ 
put plane. As the correlated areas in the crystal out¬ 
put plane are radially symmetric and practically do not 
change with intensity (see below), we can estimate the 
number of transverse modes also by the squared ratio 
determined from the appropriate widths in radial 
direction. As the curves in Fig. [5] confirm, all these quan¬ 
tities give reasonable numbers of modes of the analyzed 
twin beam close to the entanglement dimensionality Kkp. 

In the wave-vector transverse plane, decrease of entan¬ 
glement dimensionality with the increasing pump 
power Pp is explained by decrease of the width Aug^k 
of intensity profile in radial wave-vector direction (see 
Fig. [H]) accompanied by increase of widths and 

AA^ ^ of amplitude auto-correlation functions in ra¬ 
dial and azimuthal wave-vector directions, respectively 
[3 m, . The ring in the transverse wave-vector plane 
formed by the signal photons thus becomes narrower 
with the increasing pump power Pp, as confirmed by 
the radial signal-field intensity profiles Ug^k plotted in 
Fig. [9][b). The behavior of intensity profile Ug^k and in¬ 
tensity auto- (As,fc) and cross-correlation (Cg^k) functions 
in radial wave-vector direction (see Fig. ITU)) resembles 
that found in the frequency domain. Also here the auto¬ 
correlation functions Ag^k are broader than their cross¬ 
correlation counterparts Cg^k for low intensities, but they 
approach each other for more intense twin beams (see 
Fig. [TUll . This behavior follows from qualitative similar¬ 
ity of mode profiles in both variables. We remind that an 
Z-th mode in radial wave-vector direction has I maxima 
and I — 1 zeroes in its intensity profile. The behavior of 
auto- (As,(^) and cross-correlation (Cg^^) functions in the 
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FIG. 9. (a) Width An^^k of radial signal-field intensity pro¬ 
file as a function of pump power Pp and (b) radial inten¬ 
sity profile na,fc for Pp = 1 x 10“^ W (plain curve) and 
Pp = 2 X 10 “^ W (solid curve with *); Wp = lx 10 “® m, 
AAp = 1 X 10“^° m. Intensity profile ris,* is normalized such 
that f (IksTiB^kiks )/= 1- 



FIG. 10. (a) Widths of signal-field intensity auto¬ 

correlation function (solid curve) and ACs,k of intensity cross¬ 
correlation function (dashed curve) in radial wave-vector 
direction as they depend on pump power Pp. In (b), 
signal-field intensity auto-correlation function Al f.{Sks) = 
■Aa,k{ks + 5kB,k'^)/Aa^k{ks,k°) and cross-correlation function 
Cl ki^ks) = Cs,k{ka -I-(5fcs, fef’)/C's,fe(fcs, fc?) are plotted for 
Pp = 1 X 10“^ W (plain curves) and Pp — 2 y. 10“^ W 
(nearly coinciding curves with *); lOp = 1 x 10~® m, AAp = 
1 X 10”“ m. 



(a) (b) 

FIG. 11 . (a) Widths AAs,ip of signal-field intensity 

auto-correlation function (solid curve) and ACs,^ of in¬ 
tensity cross-correlation function (dashed curve) in az¬ 
imuthal wave-vector direction in dependence on pump power 
Pp. In (b), signal-field intensity auto-correlation function 
= ^s,^(</’s V5s)Ms,¥.(v5s, V^s) valid for an arbi¬ 
trary ips and cross-correlation function = Ca,Lp{f^ + 

given for an arbitrary yj? = ips + tt 
are shown for Pp = 1 x 10 ~^ W (plain curves) and Pp = 
2x 10 “^ W (nearly coinciding curves with *); uip = 1 x 10 “® m, 
AAp = 1 X 10 “^° m. 




Pp (10^" W) Fp (10“"* m) 

(a) (b) 

FIG. 12. (a) Width Ans,r of radial signal-field intensity profile 
in the crystal output plane of as a function of pump power 
Pp and (b) radial intensity profile ns,r for Pp = 1 x 10~^ W 
(plain curve) and Pp = 2 x 10“^ W (solid curve with *); 
Wp = 1 X 10~® m, AAp = 1 X 10“^'^ m. Intensity profile ns,r 
is normalized such that draVsns^rirs) = 1/2. 


azimuthal wave-vector direction is similar to that found 
in the radial wave-vector direction (see Fig. HH). 

Contrary to the transverse wave-vector plane, decrease 
of entanglement dimensionality Krip with the increas¬ 
ing pump power Pp manifests itself solely by decrease 
of the width Aug^r of radial signal-field intensity profile 
in the crystal output plane (see Fig. IT^ . Whereas the 
width Aus^r of radial signal-field intensity profile coin¬ 
cides with the width of pump beam for low-intensity twin 
beams (^ . it is narrower for more intense twin beams. 
This is explained by more intense amplification of the 
modes localized close to the pump-beam center relative 
to those occurring at the tails of the beam. Widths of 
auto- {AAg^r and AAg^pj) and cross-correlation {ACg^r 
and ACs^pj) functions as well as their shapes are nearly 
identical in the crystal output plane, as shown in Fig. 1131 
Moreover, they do not practically depend on the pump 


power Pp. 

Whereas both auto- and cross-correlation functions in 
the frequency, time and transverse wave-vector domains 
have compact shapes, long tails and oscillations are char¬ 
acteristic for the correlation functions in the crystal out¬ 
put plane (see Fig. [T5)) (^ . This stems from a different 
mode structure found in this case and discussed below. 
The analysis has shown that the correlation functions 
As^rijj and Cg^rip are rotationally symmetric and more- 
less independent on the position inside the emission disc. 
This originates in the used experimental configuration in 
which ACs^k/k^^ ~ 0.01. This value is so low that it 
does not allow to develop variations with the varying po¬ 
sition in the crystal output plane. We note that photon 
pairs emitted at the crystal end contribute to the center 
of correlation functions, whereas photon pairs generated 
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FIG. 13. 


function 


(a) [(b)] Signal-field intensity auto-correlation 
(solid curve) and cross-correlation 
function [CT,^] (dashed curve) are shown for 

Pp = 1 X 10 ”^ W (both curves nearly coincide); Al^riTs) = 
^s,7-(rs,'rg)/Ag,r(rs, Tg), C'g^j.(rg) = C's,r(rs, rf )/C's,7-(rg, rf), 
^s.v-Ws) = + 5'ips,i>°-,r°)/As,i,{'ip°,ip°-,r°), 

+ 5ij},tpf;r°,rf)/Cs,^iip°,ip?-,r°,rf); 
Tg = rf = 1.006 X 10 “® m; = ipf; Wp = 1 x 10 “^ m, 
AAp = 1 X 10 "^° m. 


at the beginning of the crystal are observed at the tails 
of the correlation functions. Thus, the width of 

radial signal-field amplitude auto-correlation function is 
sufficient for the characterization of coherence proper¬ 
ties (AAg,. = 2.297 X 10“® m). We note that the width 
AA^ ^ in the azimuthal direction depends on the dis¬ 
tance Tg from the disc center. It attains its maximal 
value (AAs^^if, = 2 tt) for rg = 0 m and then monotonously 
decreases with the increasing distance rg, in accord with 
the polar geometry. However, the presence of oscillations 
in the correlation functions shown in Fig. 1131 disqualifies 
the width AHg ^ (FWHM) as a suitable quantifier of the 
extension of field’s correlations. Width AHg ^ determined 
from the first moments of position and defined in the cap¬ 
tion to Fig. [5] has been found suitable in this case. It has 
also been used in the determination of the number 
of signal-field modes in the crystal output plane plotted 
in Fig.iKiFJ^p « 3.17iFg%). 

The oscillatory behavior of correlation functions and 
their independence on pump power Pp originate in the 
form of radial modes Us.miifs) and fti,mz(p) given by 
the transformation from the wave-vector transverse plane 
based on the Bessel functions [see Eq. (IA9I) in Ap¬ 
pendix A]. There exist two types of modes. Modes ob¬ 
tained for the azimuthal number m = 0 have their maxi¬ 
mum at r = 0 m [see Fig. I14l a)]. They are indispensable 
for describing the central part of emission disc. On the 
other hand, modes with m 7^ 0 have zero intensity for 
r = 0 m and attain their maximal values for rg^max > 0 
[for m = 1, see Fig. im b)]. The larger the azimuthal 
number m, the greater the value of rg^max- Fixing the 
azimuthal number m, all radial modes with different ra¬ 
dial numbers I have zeroes in their intensity profiles at 
the same positions. This property leads to practical inde¬ 
pendence of the correlation functions on pump power Pp. 
As the graphs in Fig. (TD indicate, the modes Us,mi with 




Tg (10 ** m) Cg (10 ® m) 

(a) (b) 

FIG. 14. Intensity profiles |Us,mi of signal-field radial modes 
in the crystal output plane for I = 0 (curve with A), 1 (*) 
and 2 (o) assuming (a) m = 0 and (b) m = 1. Modes Us,mi 
are normalized such that drgrg|Mg,mi(rs)|^ = 1; Wp = 1 x 
10"’^ m, AAp = 1 X 10“i° m. 


odd radial numbers I have small intensities compared to 
those with even radial numbers 1. So, the modes with odd 
numbers I have to be very delocalized in radial direction 
r and their influence to the properties of twin beams is 
practically negligible. This behavior has its origin in the 
shapes of modes Ug,mi in the radial wave-vector direction 
that are close to odd functions in 6k^. 


VII. CONCLUSIONS 

We have analyzed the properties of general spatio- 
spectral twin beams in the paraxial and parametric ap¬ 
proximations. Considering their spatial and spectral de¬ 
grees of freedom in their common evolution during the 
nonlinear interaction, we have investigated the properties 
of twin beams as they depend on the pump intensity. We 
have determined auto- and cross-correlation functions of 
a twin beam in the spectral and temporal domains as well 
as the transverse wave-vector and crystal output planes 
in terms of the appropriate paired Schmidt modes. We 
have mutually compared their behavior. Whereas the 
spectral and temporal coherence and the coherence in 
the transverse wave-vector plane increase with the in¬ 
creasing pump intensity, the coherence in the crystal out¬ 
put plane is practically independent on the pump inten¬ 
sity. Whereas the spectral and transverse wave-vector 
auto-correlation functions are broader than their cross¬ 
correlation counterparts for lower pump intensities, the 
opposed is true for the temporal correlation functions. 
However, auto- and cross-correlation functions approach 
each other for higher pump intensities. 

Entanglement dimensionality of a twin beam as a func¬ 
tion of the pump intensity has been determined and com¬ 
pared with the numbers of modes derived from solely the 
signal field using either its photon-number statistics or 
widths of appropriate auto-correlation functions. The 
numbers of signal-field modes have been confirmed as 
good quantifiers of entanglement dimensionality of the 
twin beam. 
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Practical independence of auto- and cross-correlation 
functions on the pump intensity in the crystal output 
plane has been explained by the special structure of 
paired modes in this plane qualitatively different from 
the common one occurring, e.g., in the spectral or tem¬ 
poral domains. Moreover, only every second paired mode 
contributes significantly to the structure of a twin beam 
for non-collinear geometries. 

We believe that this comprehensive analysis of intense 
twin beams will stimulate further experimental investi¬ 
gations of intense twin beams. Moreover, as all spatio- 
spectral modes of a twin beam are taken into account, 
the model allows for its extension to pump intensities at 
which pump depletion is observed l45j. 
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Appendix A: Properties of twin beams in the 

transverse wave-vector and crystal output planes 

In Appendix A, we define quantities in the transverse 
wave-vector space determined by spectral averaging as 
well as quantities in the crystal output plane obtained 
after temporal averaging. This averaging provides trans¬ 
verse intensity profiles as well as intensity auto- and 
cross-correlation functions. 

The signal-field intensity profile in the transverse 
wave-vector plane is obtained as follows: 

IT'S,kip{kg ,ips) — ) T’s) Ws, L)cis(^s : T))|| 

= (Al) 

q ml 

In Eq. m, symbol Oy means spectral averaging. The 
radial signal-field intensity profile Us^k is then gives by a 
cut from the intensity profile ns^kip- 

ns,kik^) = ns,kvik^,V’s = 0)- (^2) 

Averaged signal-field intensity correlations in the 
transverse wave-vector plane are described by the fol¬ 
lowing fourth-order auto-correlation function As^kip- 

As,kviki-,ips,K-^,ip'^) = {Af : A[dl{k^,ips,u}s,L) 

X dsik^, A)] A[a|(fc/, w', L) 

X A)] :)|| 

q 

(A3) 


The signal-field amplitude auto-correlation function 
^sq kip of mode q is determined as follows: 

^s,q,kipikg i^Sikg :^s) 

= (fis(fcs ; V^SI Wg, A)ds(fcs ^))\\,q 

^Y^lrnl{ki (A4) 


Radial {As^k) and azimuthal signal-field intensity 

correlation functions are derived from Eq. (IA4I) along the 
relations: 


A^^k{ki,K^) = As^kipiki-, ^3° = 0, V3'° = 0), 

As.v(V3s, (p's) = As,kv,{ki~°, Ps, ki-°, pi). (A5) 


Similarly, intensity cross-correlations in the wave- 
vector signal and idler transverse planes are quantified 
by the fourth-order correlation function Ckip'- 


Ckipik^ ,Ps,k-^ ,Pi) 

= (A/”. Ajcff (fcg , ps, LOs, L)ds{ks )V 3 s;Ws,A)] 

X A[d'l{k^,pi,uJi,L)di{k^,pi,Wi,L)] :)|| 


-E 


^ ^ ^s.ml jk^ :Ps)i\.ml{ki iPi)UmlqVmlq 


(A6) 


Radial (Cs^k) and azimuthal intensity cross¬ 

correlation functions are easily determined from the 
cross-correlation function Ckip'- 

Cs,k{k^,k^) = Cs,kp,{k^,Ps = 0, k^, p'^ = tt), 

Cs,ip{ps, Pi) = Cs,kip{k^°, Ps, k^°, Pi). (A7) 


On the other hand, properties of the twin beams at 
the crystal output plane (near field) are described by the 
two-dimensional Eourier transform of eigenmodes t writ¬ 
ten in Eq. (HU) and defined in the transverse wave-vector 
plane. This transform applied to the radially symmetric 
geometry leaves us with eigenmodes ta [xa = Xa cos(^a), 
Va = Xa sm('!/'a)], 


ts,m/(^S5 V^s) 


(^i; f/^i) 


Us,ml{Xs) exp(*m-!/'s) 

7^ ’ 

Ui^rnljxi) exp{-im'tpi) 


(AS) 


where 


^a,ml ('^a) 




(fcj-)Jm(fcj-ra), 


a = s^i 


(A9) 


and Jm stands for the Bessel function of m-th order. 

Using eigenmodes is,mi defined in Eq. (IA8I) . the aver¬ 
aged signal-field photon flux Ig^np in the crystal output 
plane is expressed as: 

ds.rv-(»'s, V's) = 2eoc{E7\rs,'tps, L,ts)E7\rs,'ips, L,t's))\\ 

= EE \is,ml{Xs,1ps)\ ^mlq- (AlO) 

q ml 
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The corresponding radial signal-field intensity profile Ig^r 
is determined as: 


IsA'^b) = =0). (All) 

The averaged signal-field intensity auto-correlation 
function Ag^rtp in the crystal output plane is obtained 
by the formula analogous to that written in Eq. (jA3l) . 

= {2eocf{Af : A[£’(“)(rs, V's, 4) 

X EA\r^, A, L, 4)] (r', As, L, Q 

xEAKAAsAAs)] :)i| 

= \^tq,rA^BAs,AAs)\^ - 

q 

In Eq. (IA12I) . the signal-field amplitude auto-correlation 
function A^ ^ characterizing mode q is determined as: 

A,q,rA^BAs,KA's) 

= 2eoc(£’i“) A, A, L, E)eA'^ (r', As, L, t'))||,, 

“ 'y 'A s.ml{'^s,'tps)is,ml{Ts,'4’sAmlq- 


The radial (Ag^r) and azimuthal {AsA signal-field in¬ 
tensity auto-correlation functions are easily derived as 
follows: 


As,r(rs,r') = Ag,r^(rg, V's = 0 ,r', = 0), 

As,^{A,As) = A,r'p,{AAs,r°A's)- (A14) 

Einally, averaged intensity cross-correlations between 
the signal and idler fields are described by the following 
fourth-order cross-correlation function: 


CrA'^B,A,n,Ai) 

= {2tocf{N : /A[EA\rB,A,LAs)AA\r^A,s,LAA\ 

X I\A[~\'ruAu LAi)EA\n,A„ LAA :)|| 


= E 


y ^ is.qAs, V^s)^i,g(r'i, 'ipi^UmlqAnlq 


ml 


(A15) 


The corresponding radial (CsA and azimuthal {CsA in¬ 
tensity cross-correlation functions are defined as: 

C's,r(rs,ri) = Cs,ri,{rsAs = 0,riAi = 0), 

AAA, A,) = Cs,rA'^° ,A,A Ai)- (Aie) 
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